1/1 


7 - 

AD  A  130  931  AIR  CODE  IMPROVEMENT  S( U I  SCIENCE  APPLICATIONS  INC  LA 
JOLLA  CA  V  A  WOOLSON  6T  AL.  01  AUG  82 
SA I  - 1 33* 02 - 232 ■ LJ  DNA  6186F  DNAOO 1  -  79  C- 0429 

UNCLASSIFIED  F/G  4/1  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  0UNCAU  OF  STANDARDS  -  i»63  -  A 


\ 

\ 


I 


l 


«A1 30931 


/t?-?3o//& 


ATR  CODE  IMPROVEMENTS 


DNA  6186F 


W.  A.  Woolson 
L.  Huszar 

Science  Applications,  Inc. 

P.  O.  Box  2351 

La  Jolla,  California  92038 


1  August  1982 


Final  Report  for  Period  15  October  1979—31  December  1979 


CONTRACT  No.  DNA  001-79-C-0429 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THIS  WORK  WAS  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODE  B325079464  V99QAXNA01015  H2590D. 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington, DC  20305 


DTIC 


l  4 


83  06 


a 


Destroy  this  report  when  it  is  no  longer 
needed.  Do  not  return  to  sender. 

PLEASE  NOTIFY  THE  DEFENSE  NUCLEAR  AGENCY, 
ATTN:  STTI ,  WASHINGTON,  D.C.  20305,  IF 
YOUR  ADDRESS  IS  INCORRECT,  IF  YOU  WISH  TO 
BE  DELETED  FROM  THE  DISTRIBUTION  LIST,  OR 
IF  THE  ADDRESSEE  IS  NO  LONGER  EMPLOYED  BY 
YOUR  ORGANIZATION. 


I 


I 

! 

i 

i 

i 

i 

UNCLASSIFIED 


security  classification  of  this  page  gMw  o mm 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  OOVT  ACCESSION  NO. 

DNA  6186F  /&>  ^ 

S.  REORIENT**  CATALOG  NUMBER 

4.  TlTL t<m4Su*Mm) 

ATR  CODE  IMPROVEMENTS 

Final  Report  for  Period 

15  Oct  79  -  31  Dec  79 

RERFORMING  ORG.  REPORT  HUNRER 

SA I- 133-82- 232-LU 

r  aotmoucj 

W.  A.  Woo Is on  and  L.  Huszar 

4  CONTRACT  or  GRANT  NUMBERS 

DNA  001-79-C-0429 

».  performing  organization  name  and  aooress 

Science  Applications,  Inc 

P.0.  Box  2351 

La  Jolla,  California  92038 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  *  WORK  UNIT  NUMWERS 

Sub task  V99QAXNA010-15 

tl.  CONTROLLING  OFFICE  NAME  ANO  AOOREM 

Director 

Defense  Nuclear  Agency 

Washington,  D.C.  20305 

It.  REPORT  DATS 

1  August  1982 

11.  NUMBER  OF  PAGES 

22  ! 

M.  MONITORING  AGENCY  name  S  AOORESSfff  AINnr t  trmm  CmmtrmUin§  Office) 

IS.  SECURITY  CLASS,  (ml  ih,t  rap mn) 

UNCUSS  IF  IED 

If*.  OECLAtSI  Fl  CATION/ DOWNGRADING 

N/fiMSTnce  UNCUSS IFIED 

IE.  OlST  RIRuTION  STATEMENT  (ml  IK  It  Km^mtt) 

Approved  for  public  release;  distribution  unlimited. 

This  work  was  sponsored  by  the  Defense  Nuclear  Agency 
under  RDT&E  RMSS  Code  8325079464  V99QAXNA01015  H2590D. 

IE.  KEY  WO  NOS  (C ANNin  «i  NrtrM  If  n«e««««ry  IPm  ttty  ly  Meet  wmNffj 

Air  Radiation  Transport  Air  Density 

Delayed  Radiation  Air  Pressure 

Atmosphere  Air  Temperature 

Fission  Isotope 

M.  AWTRACT  fCANNw  «  rntn  If  «N  ttmitty  *r  nh*  mRn; 

Modifications  were  made  to  the  ATR  code  which  will  Improve  the 
delayed  radiation  dose  model,  permit  user  specification  of  the 
air  density  profile,  and  Incorporate  additional  user  convenience 
comnands. 

00  ijSTt*  M73  «omo*»  of  » »ov tt  IS obsolete  UNCLASSIFIED 


ttcuniTv  cuSiiSclTiSi  or  Twit  page  (RMr  85?Bino 


TABLE  OF  CONTENTS 


Section  Page 

1  INTRODUCTION .  3 

2  ATR  COMMAND  MODIFICATIONS  AND  ADDITIONS .  5 

3  MODIFICATION  TO  ATR  FOR  ALTITUDE  SPECIFICATION .  7 

REFERENCES .  17 


1 


SECTION  I 


INTRODUCTION 

This  document  describes  some  modifications  to  the  ATR  code  which 
improve  the  prediction  of  delayed  radiation  dose  by  incorporating  the  varia¬ 
tion  in  fission  product  source  strength  with  the  fissioning  isotope  concen¬ 
trations,  permit  the  user  more  flexibilty  in  specifying  the  air  density,  and 
incorporate  additional  user  convenience  comnands. 

The  delayed  radiation  modification  task  was  commenced  on  the  recom¬ 
mendation  of  Jess  Marcum  (.Reference  1).  During  the  course  of  this  effort, 
whose  modest  scope  was  to  change  only  the  source  output  normalization  for 
each  fission  isotope  (^®U,  ^®U,  ^Pu)  and  to  add  comnands  to  permit 
specification  of  relative  concentration  of  each  isotope,  interest  developed 
to  advance  the  state-of-the-art  of  delayed  radiation  phenomenology  by  analyz¬ 
ing  in  detail  the  time  energy-dependent  variation  of  the  fission  product 
sources  from  each  isotope  as  well  as  the  hydrodynamic  enhancement  effects, 
fireball  rise,  and  spatial  distribution  of  sources  within  the  fireball. 
Consequently,  a  more  comprehensive  analysis  was  conducted  as  a  separate  re¬ 
search  project.  The  current  ATR  code  CATR/ISO)  with  the  improved  delayed 
radiation  module  reflecting  this  more  ambitious  scope  of  work  is  reported 
in  a  separate  document  (Reference  2). 

The  additional  user  convenience  commands  were  added  at  the  request 
of  the  U.S.  Army  Nuclear  and  Chemical  Agency  (USANCA)  in  memo  to  Dr.  David 
Auton  in  September  1979.  These  additional  comnands  are  described  in  Section 
2  of  this  report. 

The  air  density  variation  with  altitude  In  ATR  was  based  on  the 
U.S.  Standard  Atmosphere  (USSA)  and  did  not  permit  modifications  for  tempera¬ 
ture  variation  and  the  effects  of  ground  on  the  temperature  variation.  Thus 
the  air  density*  at  say,  6  Kft  above  sea  level  was  the  same  for  a  ground 
elevation  of  500  feet  or  5  Kft  above  sea  level.  In  principal,  an  ATR  user 
could  obtain  the  correct  ground  correction  factors  and  air  mass  penetration 
from  source  to  target,  but  the  procedure  was  very  cumbersome.  A  revised  air 
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density  model  was  added  with  some  new  commands  during  the  course  of  this 
effort  which  facilitates  the  detailed  specification  of  the  atmosphere  for 
ATR  calculations.  This  is  described  in  Section  3  of  this  report. 

Pending  more  substantive  modifications  to  the  ATR  code  which 
materially  effect  the  results  output  to  user,  the  changes  described  here 
will  not  be  released  to  avoid  proliferation  of  versions  of  the  code. 
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SECTION  2 


ATR  COMMAND  MODIFICATIONS  AND  ADDITIONS 

Because  of  the  wider  acceptance  by  the  community  of  tissue  kerma  as 
a  more  relevant  neutron  dose  response  than  the  Henderson  tissue  dose  response, 
the  computation  of  total  tissue  dose  which  sums  the  neutron  and  various  gamma- 
ray  components  will  now  use  tissue  kerma  for  the  neutron  component  replacing 
the  previously  used  Henderson  tissue  dose  (Reference  3). 

An  new  command  sequence  has  been  added  to  allow  the  user  a  choice 
of  ground  correction  factors.  This  sequence  is 

*C0R£  *  Turn  "off"  correction  factors,  that  is,  no  ground 
corrections  are  applied. 

*C0R1  =  Use  "old"  correction  factors.  That  is,  use  the 
correction  factors  embodied  in  ATR  Version  4  as 
described  in  Reference  3.  These  are  based  on 
calculations  by  Pace  (Reference  4). 

♦CORN^  =  Use  "new"  correction  factors.  That  is,  use  the 
correction  factors  in  ATR  Version  4.1  described 
in  Reference  5.  These  are  based  on  calculations 
by  Gritzner  and  Straker  embodied  in  TWEEDEE  Code 
(Reference  6). 

The  term  "old"  and  "new"  reflect  a  chronological  artifact  and  do  not  imply 
quality  distinctions.  The  default  command  is  CORN.  Use  of  the  above  com¬ 
mands  sets  a  switch  for  the  requested  correction  factors  which  remain  in 
effect  until  superseded  by  another  command  in  the  sequence. 

A  new  command  has  been  added  to  allow  the  user  to  specify  dose 
transmission  factors  for  individual  dose  components.  The  command  is: 

*TRAN/X/f 

where  X  *  N  for  neutrons,  NG  for  secondary  gamma  rays,  G  for  prompt  gamma- 
rays,  and  FP  for  fission  product  gamma-rays;  and  f  *  the  dose  transmission 
factor.  The  factor  is  applied  to  all  dose  calculations  and  is  operative  in 
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the  application  of  *CONTRAINT  coeaand.  The  default  value  Is  f  »  1.  The 
comand  renal  ns  In  effect  until  superseded  or  a 

•DEL  TRAN 

comand  Is  given. 


SECTION  3 


MODIFICATION  TO  ATR  FOR  ALTITUDE  SPECIFICATION 

The  objective  of  this  modification  is  to  permit  the  ATR  user  to 
specify  the  air  density  and  density  gradient  needed  by  ATR  with  greater 
flexibility.  The  basic  default  atmosphere  is  pegged  to  U.S.  Standard  Atmos¬ 
phere  (USSA)  used  in  previous  versions. 

The  air  density  is  governed  by  the  gas  equation: 

P  * 

P  * 

P  * 

R  - 
T  » 

Mo  * 

The  molecular  weight  varies  with  humidity  and  with  altitude  above 
A  molecular  temperature  T*  is  defined  so  that: 

T*  -  28.9644  jj  C2) 

which  facilitates  treatment  of  varying  molecular  weight. 

The  hydrostatic  equation  governs  the  altitude  pressure  vatiation. 

It  is: 


where 


86  km. 


o  RT 

P7T 


(1) 


2 

air  pressure  (newton/m  *  0.01  millibar) 
density  (kg/m^) 

universal  gas  constant  ^8.31432  x  10 
temperature  (°K) 

average  molecular  weight  ^28.9644 


newton  m 
kmol  ok 


■) 


dP  •  -gpdZ  (3) 

2 

where  g  ■  gravity  acceleration  g(Z-O)  ■  9.80665  m/sec 
Z  -  geometric  altitude  Cm) 
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A  geopotential  altitude  H  Is  defined  as: 


gQ  dH  »  gdZ 


(4) 


so  that  the  variation  of  gravity  can  be  eliminated.  Since  the  variation 
of  gravity  with  respect  to  Z  is: 


(5) 


where  r„ 
0 


6.356766  m,  H  is  related  to  Z  by 


H 


roz 


or 


(6) 


Before  the  gas  equation  (1)  and  the  hydrostatic  equation  (3)  can 
be  solved,  the  temperature  dependence  with  altitude  must  be  assumed.  The 
U.S.  Standard  Atmosphere  provides  a  temperature  model  linear  with  geo¬ 
potential  over  seven  regions  up  to  86  km.  Thus 

T*  -  T*  +  L1  (H  -  H.  <  H  <  (7) 


for  seven  regions  given  the  table  below. 


1 

Hi 

(km) 

L1 

(°K/km) 

o 

•r* 

1- 

P1 

(mbar) 

1 

0  and  below 

-6.5 

288.15 

1.01325+3 

2 

11 

0.0 

216.65 

2.2632+2 

3 

20 

1.0 

216.65 

5.4748+1 

4 

32 

2.8 

228.65 

8.6801+0 

5 

47 

0.0 

270.65 

1.1090+0 

6 

51 

-2.8 

270.65 

6.6938-1 

7 

71 

-2.0 

214.65 

3.9564-2 

84.5 

187.65 

3.9814-2 
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Equations  (9)  and  (10)  define  the  U.S.  Standard  Atmosphere. 

The  P.  (pressures  at  the  seven  boundary  layers)  are  calculated  from  the 
bottom  to  the  top  by  setting  *  1.01325(.+3)  mbar  and  calculating  Pj+^ 
for  H^+j  in  either  (9)  and  (10).  This  treatment  is  valid  below  85  km. 

The  default  atmosphere  in  ATR  is  the  U.S.  Standard  Atmosphere, 
but  the  user  can  now  modify  It  with  the  new  command  structure  by  speci¬ 
fying  other  conditions  at  a  single  altitude,  or  the  user  can  force  a  den¬ 
sity  profile  on  a  single  altitude  band  within  this  standard  atmosphere  by 
inputting  pairs  of  altitude  and  density  values. 
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The  single  point  altitude  change  commands  are: 

*ELEV/UNIT/HP  -  specification  elevation  -  single  point  where  user  will 
specify  modifications  to  the  USSA.  Default  value  of  HP 
is  ground  elevation.  UNIT  =  M,  FT,  KM,  KFT,  MI, 
default  is  meters. 

*TEMP/UNIT/TP  -  temperature  at  HP.  Default  is  calculated  from  input  PP 
and  DP  or  15°C.  Unit  *  C  (°C),  K  (°K),  or  F  (°K), 
default  is  °C. 

*PRES/UNIT/PP  -  pressure  at  HP.  Default  is  calculated  from  input  TP  and  PP 
or  USSA  value  at  HP.  UNIT  =  MB  (millibar),  TR  (torr), 

HG  (in  Mercury),  PSI  (lbs/in2),  LBF  (lbs/ft3), 
default  is  MB. 

* DENSE/UN IT/DP  -  density  at  HP.  Default  is  calculated  from  TP  and  PP. 

UNIT  *  KGM  (.Kg/m3),  GMC  (gm/cm2),  LBF  (lbs/ft3) 
default  is  GMC. 

NOTE:  An  input  DP  and  TP  will  over  ride  an  input  PP. 

The  single  point  modification  of  USSA  is  performed  in  the  follow¬ 
ing  way. 

Equation  (1)  defines  a  triplet  between  TP,  PP,  and  DP  at 
altitude  HP.  If  all  three  are  specified,  only  TP  and  DP  are  used  and 
the  input  value  of  PP  is  ignored.  If  any  two  are  input,  the  third  is 
calculated  from  equation  (1).  If  none  or  only  one  are  specified,  the 
default  values  are 

TP  =  15°C 

PP  =  USSA  value  from  equations  (.9)  or  (10) 

DP  ■  calculated  from  equation  (1)  and  TP,  PP  values. 

This  defines  the  atmosphere  only  at  HP. 
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Define  the  geopotential  altitude  HG  by: 


HG 


r  HP 
o 

r  +  HP  • 
o 


(11) 


and  to  calculate  the  rest  of  the  atmosphere,  if  needed,  the  values  of  and 
P..  (Table  1)  nearest  to  the  altitude  nearest  HG  are  adjusted.  This  requires 
adjusting  two  L. 's  so  that  the  temperature  dependence  is  continuous.  In  the 
lowest  altitude  band  we  only  change  L^.  The  procedure  for  modifying  T^,  L. , 
and  Lj+1  is  illustrated  in  Figure  1  for  two  examples.  The  figure  presents 
the  USSA  temperature  vs.  geopotential  altitude. 

First  consider  a  user  specifying  a  point  Pj  in  Figure  1  correspond¬ 
ing  to  values  HG,  TP,  PP,  and  DP.  HG  is  closer  to  H^  than  H2.  Therefore 
Tj,  P^,  and  will  be  modified  by  linear  interpolation 


T 


H2  -  HG 

npnq 


P£  must  be  calculated  from  equation  (9)  (use  equation  (10)  if  ■  0),  via 


P1  3  PP  [jP  +  L'Thj  -  HG)  j 

Now  with  Pj,  T. ,  and  Lj  redefined  to  P£,  Tj,  and  LJ  equation  (9)  or  (10)  can 
be  used  to  get  any  density. 

In  the  second  case  In  Figure  1  the  user  specified  altitude  which 
is  not  near  the  ground.  The  procedure  Is  similar  except  that  T^,  P^,  Lj, 
and  must  be  modified  in  order  to  assure  temperature  continuity.  Here 
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HG  Geopotential  Altitude  (km) 


T  oc 


Figure  1.  USSA  temperature  vs.  geopotential  (solid  line). 
Two  examples  of  user  modifications  shown  as 
dotted  lines. 


(TP  -  T3) 


Now  we  can  specify  the  general  case  procedure.  First  we  calculate 
HG  x  and  test  HG  against  H.  to  get  1.  Then  we  test  whether  HG 

Tq  +  Hr  1 

is  above  or  below  the  midpoint  of  that  1: 

If  HG  <  (Hi  +  Hi+1)/2,  we  modify  T^ ,  Pi#  L^,  and  L^; 

If  HG  >  (H1  +  H1+1)/2,  we  modify  T1+1,  P,+1«  and  L1+r 


Now,  if  HG  <  (Hi  +  Hi+1)/2,  we  have 

H..,  -  HG 

T1  *  T1+l  "  Hi+1  -  fl1  (T1+1 


TP) 
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and 


l; 


t;  -  Vi 


f-l  *  H,  -  Hf_ 


(not  needed  if  i  =  1) 


and 


1 


However,  if  H6  >  (Hi  +  H.+1)/2,  we  have 

HT  -  Tj 

TT  ,  -  Tj  +  u - r  (TP-^) 

i+1  1  H.+1  -  Hi  1 


and 


and 


1  Hi+1  Hi 


,  .  Ti+2  "  Ti+i_ 

1+1  H1+2  -  H1+1 


(not  needed  if  1  *  6) 


and 


pp 


g  M. 
ao  o 

TT7 


r _ s - iv 

[TP  ♦  l,'  (H)+1  -  HS)  J 
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If  the  user  wishes  to  over  ride  the  density  profile,  the  comnands 

are: 

ALT/UNIT/HRj  HRg....  -  array  of  altitudes  to  specify  a  density  profile,  no 

default.  Units  same  as  ELEV  command. 

HR.  <  HRi+1. 

RHO/UNIT/HDj  HDg....  -  array  of  densities  corresponding  to  HR,  no  default. 

Units  are  same  as  DENSE  command. 

When  the  ATR  user  inputs  pairs  of  (HR^ ,  DR^)  to  specify  his  own 
density  profile,  these  override  the  previous  atmosphere  in  the  range  in¬ 
cluded  in  HRj  to  HRN.  Interpolation  between  HR^  and  HR^+^  assumes  an 
exponential  atmosphere  and  is  performed  by 


P(Z) 


However  if 


pi+l  "  P1 
pi 


<  0.001  linear  Interpolation  is  used  via 


p(Z) 


n  ,  p1+l  ~  P1 

P1  +  (Hki+1  -  HR1 ) 


(Z  -  HRi ) 
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